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Fluid forces on. a rotating centrifugal impeller in whirling motion were 
studied experimentally. A two-dimensional impeller installed in a parallel- 
walled vaneless and vaned diffuser whirled on a circular orbit with various 
positive and negative angular velocities. The results showed that the fluid 
forces exert a damping effect on the rotor in most operating conditions, but 
become excitatory when the impeller operates at very low partial discharge while 
rotating far faster than the whirl speed. 

The fluid forces were also expressed in terms of mass, damping and stiff- 
ness matrices. Calculations were conducted for impellers with the same geometry 
and whirl condition as those in the experiment. Quantitative agreement was ob- 
tained especially in positive whirl. 

INTRODUCTION 

As to fluid forces acting on a centrifugal impeller of pumps and compres- 
sors, many studies have been done on the steady lateral forces (radial thrust), 
which occur either by azimuthal ly asymmetric configuration of volute casing or 
by misalignment of impeller and casing center. To the contrary, limited contri- 
butions have been made to the unsteady part of lateral forces, which is induced 
by the whirling motion of rotating impellers. This is the very information 
sought for the precise analysis of rotor stability and seismic response. 

Experimental study on this subject is quite few* Two projects are pre- 
sently under way as far as known to the authors, one at California Institute of 
Technology and the other at University of Tokyo. 

The latter, i.e. the present experimental study is the direct counterpart 
to the theoretical calculation reported at the 1st Workshop [ref. 1]. Therefore 
the geometries of the test impeller and casing were determined to be similar to 
those of sample calculations. This is the reason why this experiment adopted 
deliberately much simpler configurations than those of actual pumps and compres- 
sors. The project at CIT seems to lay more stress on the practical side than we 
intend, because the adopted impeller/casing configuration is similar to those of 
pumps for industrial application. The progress of CIT project was reported 
twice at the 1st and 2nd Workshop by Acosta et al. [ref. 2, 3]. 


* This study was supported by Grant in Aid for Scientific Research by the 
Ministry of Education, Technical Research Institute for Integrity of Structures 
at Elevated Service Temperatures and Hitachi, Ltd. Messrs. S. Yanagisawa, 

K. Tomita, J. Hanawa, K. Kawakami and C. Kato contributed to this study as their 
Master theses. 
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TEST FACILITY 


The test pump is a single-stage, vertical shaft, centrifugal type. The 
mechanism of the forced whirling motion is illustrated schematically in fig, 1. 
The impeller is supported overhung by two self -alignment roller bearings, which 
are mounted separately on the upper and lower whirl plates. Fluid forces acting 
on the impeller are evaluated from the reaction force to the lower bearing. 

As seen from fig. 1, the outer race of the lower bearing is attached to the 
lower whirl plate by four load cells arranged with 90 degree angle difference 
(x- and y-direction) . 

Fig. 2 shows the relative arrangement of impeller and diffuser. In order 
to avoid the generation of steady lateral forces, the flow passage was designed 
to ensure axisymmetric flow by means of resistance screen at the diffuser exit 
and two symmetrically arranged discharge pipes for instance. Test was conducted 
mainly for vaneless diffuser. Additional test was also done for vaned diffuser 
with 2, 4 and 8 guide vanes, though the corresponding theory had not been estab- 
lished for comparison. Two impellers, A and B, with different vane angle were 
tested. The rotational speed of the impeller shaft was kept constant (8.8 Hz), 
while the whirl speed was varied widely from negative to positive whirl. Both 
pump and eccentric drive shafts were equipped with rotary encoders to detect 
phase angle. Principal specifications of test pump are listed in table 1. 


DATA EVALUATION 


Calibration 

The load cells were calibrated dynamically by a rotating disc with known 
unbalance mass. The conversion factors (a = impeller force (N)/load cell output 
(V) ) were determined for the x- and y-direction separately. The linearity of 
the cells was satisfactory and all cross terms, output of y-cells to x-directed 
force for instance, remained less than two percent of their principal terms. 


Data Aguisition 

Outputs of load cells were A/D converted and stored by sampling signals 
generated at every 10 degree on whirling orbit. For experiment to check the 
statistical character of the data, ensemble average and rms were evaluated from 
512 sampled data at the same condition. In other cases 200 to 300 data were 
collected to give an ensemble average with intended confidence interval. In the 
experiment with vaned diffuser, samplings were made when the phase angle of 
whirl and impeller shaft, 0^ and 0^ (cf . fig, 4(a)), coincided simulta .eously 
with the preset values. 

Evaluation of Fluid Forces on Impeller Vanes 

Since the theoretical results which are to be compared with the experiment, 
consider only fluid forces acting on the impeller vanes, and exclude all other 
forces on shroud, shaft, seal etc, it is necessary also in the experiment to ex- 
tract pure vane force from the resultant reaction force measured at the lower 
bearing. For this purpose, four experiments were carried as shown in fig. 3, 
that is, experiments with whirling impeller in water and in air, and with whirl- 
ing dummy rotor in water and in air. 

The fluid force on the impeller vanes, F(F X , F ), C an be evaluated by the 
relation, ^ 


F = a 




+ E. ) 


( 1 ) 


where E. (E , 
1 x 


E^) denotes load cell outputs in the above mentioned four sub- 
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experiments. The following interpretation will also help understand the meaning 
of eq. (1): 

a(£ 1 ~E 2 ^ ' fluid forces on whirling components submarged in water. 

a(E^“E 4 ) * fluid forces on whirling components submarged in water except 
vanes, i.e. those on shrouds, shaft and seal (see fig. 1). 

The reaction of driving shaft torque on the lower bearing was calculated 
and the influence was compensated in the final evaluation of impeller force. 

As seen from fig. 4 (a), fluid forces on whirling impeller are dependent of 
rotational speed u> , whirl speed ft , eccentricity € , phase angle of whirl 6^ 
and rotation 0 and flow rate Q. The dependence on 0 disappers in vaneless 
diffuser. 

Test Condition 

The eccentricity was kept constant (£ = 1.5 mm) , since linearity could be 
expected for such small displacement. From the rotor dynamic consideration the 
whirl speed ft corresponds to the first damped natural frequency in bending mode 
of the rotor (critical speed). Although most troubles associated with non-syn- 
chronous rotor vibration occur at supercritical rotational speed, ft/UK 1, the 
test was done at whirl speed ratio ft/co up to 1.3 both for positive and negative 
whirl. Because impeller forces are sensitive to the operating condition of the 
impeller, discharge rate was varied from shutoff to maximum. 

Presentation of Fluid Forces 

Fluid forces F acting on the impeller vanes are divided into two compo- 
nents as illustrated in fig. 4 (b) . The component normal to the orbit and di- 
rected outward is called radial force F , while the component parallel and op- 
posite to the clockwise (positive) whirl is called tangential force Fg. 

A negative sign means the direction is opposite to the above definition. 

TEST RESULTS AND DISCUSSION 


Pump Characteristics 

The hydraulic performance of the test pump (impeller A and B) at steady op- 
eration is plotted in fig. 5 in non-dimensional form, in which flow coefficient 
<1> and head coefficient ^ are defined by (see notation in table 1) 

4) = Q/2TTr 2 bu 2 • u 2 = lor 2 (2) 

\J) = H / (u 2 2 /2g) 

Since impeller B has smaller exit vane angle than that of A, it has larger work 
input and hence higher head. Guide vanes have logarithmic spiral form and the 
vane angle is set to have no loading at design (shock-free) discharge of the im- 
peller. Because of this design principle, the influence of guide vanes on the 
steady performance is not remarkable. Design flow coefficient, at which the 
relative inlet flow to the impeller vanes becomes shock-free, is denoted by ■<!> , 

g The machine Reynolds number of the test pump. Re = u^d^/^# WaS aix)ut 3.4 S x 
10 in all test cases. 


Character of Unsteady Fluid Forces 

Fluid forces acting on whirling vanes, F , are evaluated from four forces 
obtained by four sub-experiments (see eq. (1)). The relative magnitude and di- 
rection of these four forces are illustrated in fig. 6 for a typical test con- 
dition (impeller A, , £ ~2 mm, ft/U)=0.83 and 0^=0) . and F^ measured in air 

are caused by inertial centrifugal force and thus have only radial component. 
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The effort to reduce the mass of whirling components by using aluminum helped 
keep the level of fluid force substantially larger than that of inertia force. 

Outputs of load cells oscillate primarily with the whirl speed Q, but con- 
tain also fluctuations caused by mass unbalance of impeller, rolling noise of 
bearings, noise of measuring system, turbulence of flow and so on. From the 
preliminary studies, it became clear that the largest and predominant source of 
data scattering resulted from the fluctuation of flow field in impeller and dif- 
fuser. Fig. 7 illustrates Lissajous figures of overall impeller force., when 
impeller A rotates in water without whirling motion. As seen from fig. 7 (a), 
fluctuation of measured force is rather trivial at design, i.e. shock-free dis- 
charge condition. As the discharge decreases, the fluctuation gets larger and 
at shutoff the fluctuation becomes quite violent as seen from fig. 7 (b) . 

The reason is obviously the separation and large scale turbulence of the flow in 
impeller. The fluctuation of measured force can therefore be attributed mostly 
to the inherent flow turbulence. 

Fig. 8 shows fluid forces on whirling impeller vanes at 9 locations of the 
orbit (impeller A, e-1.5 mm , vaneless dif.) for the combination of two flow 
rates (design and shutoff) and two whirl speed ratios (positive and negative). 
Each vector indicates ensemble average of 512 sampled data, while the magnitude 
of fluctuation is indicated by a circle with rms as its radius. Fluid forces 
are nearly axisymmetric and uniform on the orbit as they should be. It is worth 
to note, that especially near shutoff the fluctuation in rms can be severalfold 
as high as its average. 

Fluid Forces on Whirling Impeller in Vaneless Diffuser 

In spite of precautions to establish azimuthally symmetric flow, the en- 
semble averaged force vectors on the orbit were not completely uniform as seen 
from fig. 8. Since the influence of fluid forces on the dynamics of rotor is 
accumulative, the force components, F and F0, are hereafter averaged over one 
cycle of whirl angle 0^ and shaft angie 0 2# and are denoted by and TFg.. These 
forces are further normalized by the following definitions? 

f r = F r /MEU) 2 , f 0 = F 0 /M£OJ 2 (3) 

2 

M = pirr 2 b 

where p is the density of the working fluid and M is the mass of fluid displaced 
by the silhouette volume of impeller vanes. 

Fig. 9 shows the dependence of .f and f q on the whirl speed ratio Q/U) from 
-1.3 to +1.3, when impeller A (3=68°) whirls with eccentricity £ =1.5 mm in vane- 
less diffuser and operates at five different flow rates from shutoff to maximum. 
The broken line in the figure indicates calculated result for shock-free entry 
condition ((N^jJ. Fig. 10 shows the corresponding result of impeller B (3=60°). 

As to radial component followings can be observed: 1) The apparent mass of 
impellers A and B is about 0.9 and 0.7 times displaced fluid mass M, respective- 
ly, when it merely oscillates without rotation. Therefore, if there is no in- 
fluence of pumping action on^he fluid forces, f r must be equal to 0.9x(R/0); 
for impeller A and 0.7x(Q/W) for impeller B. The calculated f r at shock-free 
entry agrees with this simple prediction quantitatively in positive whirl, but 
becomes almost twice as large as the prediction in negative whirl. It suggests 
that the pumping impeller in negative whirl has a certain flow mechanism to ge- 
nerate large radial force even in potential theory. 2) At design discharge, 
that is, shock-free entry condition, the measured radial forces agree with the 
theory well in tendency. In positive whirl quantitative agreement is fairly 
good, while in negative whirl the measurement is remarkably larger than the 
theory. 3) Radial force decreases as discharge decreases in positive whirl. 
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In negative whirl the tendency is quite opposite. 4) At and near Q/a)=1.0 f where 
the relative flow in impeller becomes stationary, radial force becomes negative 
at very low discharge. 5) At f2/o>=0 (stationary displacement), rotating impeller 
produces positive radial force, that is, force to the direction of impeller dis- 
placement. This means that radial fluid force has negative-spring characteris- 
tics. 6) Impeller A has slightly larger radial force than that of impeller B. 

As to tangential component: 1) Except a few test conditions, tangential 
forces are positive in positive whirl and negative in negative whirl. This indi- 
cates that fluid forces exert a damping effect on whirling impeller in most 
cases. 2) Negative tangential force emerges in positive whirl when the pump 
operates at very low discharge while rotating more than twice as fast as whirl 
speed (ft/(0<0.5). In this condition the fluid force feeds energy to the whirling 
system and can cause a self-excited oscillation, a sort of impeller whip. 

3) The agreement with theory is good in positive whirl. In negative whirl the 
measurement is again remarkably higher than the theory predicts. 4) In nega- 
tive whirl larger damping can be expected than in positive whirl. 5) The influ- 
ence of discharge rate on tangential component is generally small. 

Influence of Guide Vanes on Rotor Stability 

Since the radial clearance between impeller and guide vanes was chosen re- 
latively large (11% of radius), the hydrodynamic interference of guide vanes to 
impeller vanes was so limited that the dependence of impeller force on the phase 
angle of impeller shaft, 0 , was undetectably small. Consequently the instal- 
lation of guide vanes resulted in little difference of the mean components of 
fluid forces from those given in fig. 9 and 10. 

Guide vanes augmented the damping effect on the rotor stability to some ex- 
tent. For the convenience of stability judgement, a new normalized parameter of 
tangential force is introduced as 

f Q ' ■= Fg/MeQo) (4) 

According to this definition small f a at small whirl speed ratio ft/U) can be en- 
larged by dividing by £2. fJ has also the feature that it keeps positive 
throughout positive and negative whirl range, as far as the fluid force exerts 
damping effect* Fig. 11 shows fg* of impeller A at various flow rates, when it 
whirls in vaneless (no GV) and vaned diffuser with 2 and 8 guide vanes. In the 
case of vaneless diffuser, theoretical prediction for shock-free entry condition 
is also plotted by broken lines. From these figures it is evident that guide 
vanes reduce the range of whirl speed ratio with negative fg r , thus resulting in 
increased damping on the whirling rotor. 

Comparison with Theory by Vector Diagram 

Fluid forces at shock-free entry condition are compared with the theoret- 
ical result [ref. 1] using vector diagrams as shown in fig. 12. Normalized 

force vectors are plotted with whirl speed ratio ft/U) as their parameter on the 
curves. Theoretical and experimental curves coincide fairly well as a shape of 
curves but there is a considerable difference of the location of parameters. 

This difference could be caused either by the viscous effect neglected in this 
theory or by unaware factors hidden in the test setup and instrumentation. 

CONVERSION TO MATRIX ELEMENTS 

In the analysis of rotor dynamics it is convenient to express fluid forces 
in terms which are proportional to the displacement, velocity and acceleration 
of impeller center. As illustrated in fig. 4 (a), the location of impeller 
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center 0(x f y) is expressed by the absolute coordinate system with its origin 
at the whirl center 0*. Denoting x- and y- components of fluid forces by F and 
F , they are expressed in the following linear form; 
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where x = € cos©^, y - £ sin0^ and 0^ = -Qt. 

Because the measured fluid forces are non-linear with respect to location, 
velocity and acceleration of impeller center, it is impossible to determine the 
above matrix elements in an exact meaning. However, if the fluid forces are 
axisymmetric on the whirling orbit (independent of 0^ ) , and can be approximated 
by the equations 


F ( £ , ) = E ( F + - F fi+ F ft 2 ) 
r rO rl r2 




(6) 


Fg(e, ft ) = e ( Fg 0 + F 01 Q+ f 02 ft 2 ) 
the matrix elements can be related to the above terms as; 
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The fitting of measured data to eq. (6) can be performed by applying the method 
of least squares. 

Fig. 13 shows the fitted curve of impeller A at shock-free condition. The 
fitting was made for all data in positive and negative whirl, — 1 . 3 < f2/co <+1.3. 
Since the real fluid forces are non-linear, the fitted curve cannot represent 
the measured data satisfactorily in the whole whirl speed ratio. Especially the 
delicate change of fluid forces near fi/U) = 1 is neglected in the fitting com- 
pletely. The quality of fitting deteriorates gradually as flow rate decreases 
and the flow in impeller deviates from shock-free condition. 

Fig. 14 illustrates the fitted curves of impeller A in vaneless diffuser at 
various flow rates. As seen from the figure tangential force fg changes almost 
linearly with whirl speed ratio. This indicates that cross mass term m =-m 
is insignificant. The feature of radial force f consists in the fact X ?hat y rt 
is parabolic in general but the minimum takes place around Q/U)=0.4. This result 
infers that cross damping term b =-b plays an important role. 

In table 2 and 3 matrix elem^nts^of impeller A and B in vaneless diffuser 
are listed for various flow rates. All values are normalized by the quantities 
given in the tables. 


CONCLUSION 

Fluid forces acting on two-dimensional centrifugal impeller whirling in 
vaneless and vaned diffuser were measured and compared with the corresponding 
theoretical results. Principal findings from two tested impellers are: 

1) Fluid forces have tangential components to damp whirling motion of the 
rotating shaft in most operating condition. 
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At low flow rate near shutoff and low whirl speed ratio, i.e. at high 
supercritical rotational speed, there are possibilities that the fluid 
forces do exert negative damping to the rotor in positive whirl, a sort 
of impeller whip. 

3) Theoretical calculation can predict fluid forces qualitatively. Quan- 
titative agreement can be obtained in positive whirl. 

4) Fluid forces are expressed as elements of mass, damping and stiffness 
matrices by assuming linear relation. 

This study will be extended further to the cases in which impeller whirls 
in a volute casing. The 2nd generation test rig with improved whirling mecha- 
nism, instrumentation and versatile applicability is under construction. 
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Table 1 Principal specifications 


Impeller 

Two-dimensional, closed type 

outer diameter d^=2r 2 
inner diameter 

350 mm 

175 mm 

vane width b 

3 5 mm 

suction diameter 

160 mm 

number of vanes 

6 

vane angle to radius 3 

68° for impeller A 

(logarithmic spiral) 

60° for impeller B 

mass 

4.1 kg 

Vaneless diffuser 

Parallel walled 

exit diameter 

700 mm 

screen 

40 mesh per inch, 51% opening 

width 

7 0 mm 

Vaned diffuser 

Logarithmic spiral, 5 mm thickness 

inner diameter 

390 mm 

(20 mm radial clearance) 

outer diameter 

600 mm 

width 

60 mm 

number of vanes 

2, 4 and 8 

vane angle to radius 

75.5° 

Test condition 


rotational speed W 

525 rpm 

whirl speed ^ 

-683 to +683 rpm 

eccentricity £ 

variable^ mainly 1.5 mm 

2.25 m^/min for impeller A 
3.21 m /min for impeller B 

shock-free flow rate 
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Table 2 Matrix elements of impeller A 


matrix element 

m xx* m yy 

m .-m 
xy J yx 

b , b 
xx yy 

b ,-b 
xy yx 

k , k 
xx’ yy 

k ,-fc 
xy yx 

normalized by 

irpT2 2 b 

irpr2 2 bai 

irpr 2 2 bui 2 

$ = 0 

1.64 

0.27 

1.14 

-2.02 

-0.55 

-0.13 

°- 2 ♦,£ 

1.72 

0.35 

1.10 

-2.01 

-0.55 

-0.17 

°.6 * sf 

1.87 

0.32 

1.08 

-1.88 

-0.56 

-0.10 

♦sf 

1.86 

0.27 

0.98 

-1.33 

-0.42 

0.09 

1,3 ^sf 

1.82 

0.23 

1.09 

-1.01 

-0.43 

0.08 


Table 3 Matrix elements of impeller B 


matrix element 

m yy 

X 

a* 

1 

* 

a* 

b , b 
xx’ yy 

b ,-b 
xy' yx 

k , k 
xx yy 

k , -k 
xy' yx 

normalized by 

irpr 2 2 b 

irpr 2 2 bw 

Trpr 2 2 bta 2 

<t> = 0 

1.30 

0.21 

0.98 

-1.98 

-0.58 

-0.12 

0.2 * s£ 

1.32 

0.33 

0.94 

-1.96 

-0.57 

-0.17 

0.6 * sf 

1.26 

0.20 

1.06 

-1.80 

-0.43 

-0.13 

♦s£ 

1.34 

0.19 

1.01 

-1.62 

-0.39 

-0.05 

1.3 * s£ 

1.32 

0.05 

1.00 

-1.35 

-0.35 

0.01 
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Fig. 2 Configuration of impeller 
and diffuser 


Fig. 1 Whirling 



Experiment 1 
(Impeller in water) 


mechanism 
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Experiment 2 
(Impeller in air) 



Experiment 3 
(Dummy in water) 

Fig. 3 Four 



Experiment 4 
(Dummy in air) 
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Fig. 4 Notations 
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(a) Tangential component 





















